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The PHENIX Collaboration at the Relativistic Heavy Ion Collider has measured open heavy flavor 
production in minimum bias Au+Au collisions at = 200 GeV via the yields of electrons from 

semileptonic decays of charm and bottom hadrons. Previous heavy flavor electron measurements 
indicated substantial modification in the momentum distribution of the parent heavy quarks due 
to the quark-gluon plasma created in these collisions. For the first time, using the PHENIX silicon 
vertex detector to measure precision displaced tracking, the relative contributions from charm and 
bottom hadrons to these electrons as a function of transverse momentum are measured in Au-|-Au 
collisions. We compare the fraction of electrons from bottom hadrons to previously published 
results extracted from electron-hadron correlations in pj-p collisions at = 200 GeV and find 

the fractions to be similar within the large uncertainties on both measurements for pt > 4 GeV/c. 

We use the bottom electron fractions in Au-fAu and p-\-p along with the previously measured heavy 
flavor electron Raa to calculate the Raa for electrons from charm and bottom hadron decays 
separately. We find that electrons from bottom hadron decays are less suppressed than those from 
charm for the region 3 < pt < 4 GeV/c. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

High-energy heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Col¬ 
lider (LHC) create matter that is well described as an 
equilibrated system with initial temperatures in excess 
of 340-420 MeV [THS]. In this regime, the matter is un¬ 
derstood to be a quark-gluon plasma (QGP) with bound 
hadronic states no longer in existence as the temperatures 
far exceed the transition temperature of approximately 
155 MeV calculated by lattice quantum chromodynamics 
(QCD) |2]. This QGP follows hydrodynamical flow be¬ 
havior with extremely small dissipation, characterized by 
the shear viscosity to entropy density ratio p/s « l/47r 
and is thus termed a near-perfect fluid [I1I7H9]. 

Charm and bottom quarks (m^ ~ 1.3 GeV/c^ and 
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mb « 4.2 GeV/c^) are too heavy to be significantly pro¬ 
duced via the interaction of thermal particles in the QGP. 
Thus the dominant production mechanism is via hard 
interactions between partons in the incoming nuclei, i.e. 
interactions that involve large momentum transfer, q^. 
Once produced, these heavy quarks are not destroyed by 
the strong interaction and thus propagate through the 
QGP and eventually emerge in heavy flavor hadrons, for 
example D and B mesons. 

Early measurement of heavy flavor electrons from the 
PHENIX Collaboration in Au-I-Au collisions at RHIC in¬ 
dicated that although the total heavy flavor production 
scales with the number of binary collisions within un¬ 
certainties nniiii], the momentum distribution of these 
heavy quarks is significantly modified when compared 
with that in p-\-p collisions da [13]. These results indi¬ 
cate a large suppression for high-pr > 5 GeV/c electrons 
and a substantial elliptic flow for px = 0.3-3.0 GeV/c 
electrons from heavy quark decays. Here, and through¬ 
out the paper, we use “electrons” to refer to both elec¬ 
trons and positrons. The suppression of the charm quark 
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has since been confirmed through the direct reconstruc¬ 
tion of D mesons by the STAR Collaboration [14]. In 
Pb-I-Pb collisions at the LHC at = 2.76 TeV, simi¬ 

lar momentum distribution modifications of heavy flavor 
electrons and D mesons have been measured [HI [Hj. 
Recently, the CMS experiment has reported first mea¬ 
surements of B —)■ J/ip [T7] and b-jets [T5j in Pb-I-Pb 
collisions. In contrast to this suppression pattern found 
in Au-I-Au collisions, d-l-Au and peripheral Cu-fCu col¬ 
lisions at = 200 GeV exhibit an enhancement at 

intermediate electron in the heavy flavor electron 
spectrum [MHO] that must be understood in terms of 
a mechanism that enhances the px spectrum, e.g. the 
Cronin effect [2T|- That mechanism potentially moder¬ 
ates the large suppression observed in Au-|-Au collisions 

'J^NN = 200 GeV. It is notable that in central Au-I-Au 
collisions at y/s~^ = 62 GeV an enhancement is also ob¬ 
served at intermediate px [22] . 

The possibility that charm quarks follow the QGP 
flow was postulated early on [H], and more detailed 
Langevin-type calculations with drag and diffusion of 
these heavy quarks yield a reasonable description of the 
electron data pHt29j . Many of these theory calculations 
incorporate radiative and collisional energy loss of the 
heavy quarks in the QGP that are particularly impor¬ 
tant at high-pT, where QGP flow effects are expected to 
be sub-dominant. The large suppression of heavy flavor 
electrons extending up to pt ~ 9 GeV/c has been a par¬ 
ticular challenge to understand theoretically, in part due 
to an expected suppression of radiation in the direction 
of the heavy quarks propagation - often referred to as 
the “dead-cone” effect [5^ . 

This observation of the high-p^ suppression [SI1I32] is 
all the more striking because perturbative QGD (pQGD) 
calculations indicate a substantial contribution from bot¬ 
tom quark decays for px > 5 GeV/c [53]. First measure¬ 
ments in p-l-p collisions at 200 GeV via electron-hadron 
correlations confirm this expected bottom contribution to 
the electrons that increases as a function of px [331 ES] . 
To date, there are no direct measurements at RHIC of 
the contribution of bottom quarks in Au+Au collisions. 

For the specific purpose of separating the contribu¬ 
tions of charm and bottom quarks at midrapidity, the 
PHENIX Collaboration has added micro-vertexing capa¬ 
bilities in the form of a silicon vertex tracker (VTX). The 
different lifetimes and kinematics for charm and bottom 
hadrons decaying to electrons enables separation of their 
contributions with measurements of displaced tracks (i.e. 
the decay electron not pointing back to the collision ver¬ 
tex). In this paper, we report on first results of separated 
charm and bottom yields via single electrons in minimum 
bias (MB) Au-I-Au collisions at = 200 GeV. 


II. PHENIX DETECTOR 

As detailed in Ref. [35], the PHENIX detector was 
originally designed with precision charged particle recon¬ 


struction combined with excellent electron identification. 
In 2011, the VTX was installed thus enabling micro- 
vertexing capabilities. The dataset utilized in this anal¬ 
ysis comprises Au-I-Au collisions at ,/s^ = 200 GeV. 

A. Global detectors and MB trigger 

A set of global event-characterization detectors are uti¬ 
lized to select Au-l-Au events and eliminate background 
contributions. Two beam-beam counters (BBG) covering 
pseudorapidity 3.0 < |ry| < 3.9 and full azimuth are lo¬ 
cated at ± 1.44 meters along the beam axis and relative 
to the nominal beam-beam collision point. Each of the 
BBCs comprises 64 Cerenkov counters. 

Based on the coincidence of the BBCs, Au-l-Au colli¬ 
sions are selected via an online MB trigger, which requires 
at least two counters on each side of the BBC to fire. The 
MB sample covers 96 ± 3% of the total inelastic Au-I-Au 
cross section as determined by comparison with Monte 
Carlo Glauber models m- The BBC detectors also en¬ 
able a selection on the z-vertex position of the collision as 
determined by the time-of-flight difference between hits 
in the two sets of BBC counters. The z-vertex resolu¬ 
tion of the BBC is approximately = 0.6 cm in central 
Au-I-Au collisions. A selection within approximately ±12 
cm of the nominal detector center was implemented and 
~ 85% of all Au±Au collisions within that selection were 
recorded by the PHENIX high-bandwidth data acquisi¬ 
tion system. 

B. The central arms 

Electrons (e'*' and e~) are reconstructed using two cen¬ 
tral spectrometer arms as shown in Fig. [Ra), each of 
which covers the pseudorapidity range \r]\ < 0.35 and 
with azimuthal angle Acj) = 7r/2. The detector config¬ 
uration of the central arms is the same as in previous 
PHENIX Collaboration heavy flavor electron publica¬ 
tions [HdS]. Charged particle tracks are reconstructed 
outside of an axial magnetic field using layers of drift 
chamber (DC) and multi-wire proportional pad cham¬ 
bers (PC). The momentum resolution is (Jp/p ~ 0.7% 
0 0.9% p (GeV/c). Eor central arm charged particle 
reconstructions the trajectory is only measured for ra¬ 
dial positions r > 2.02 meters, and the momentum vec¬ 
tor is calculated by assuming the track originates at the 
Au±Au collision point determined by the BBC detectors 
and assuming 0 radial distance. 

Electron identification is performed by hits in a ring 
imaging Cerenkov detector (RICH) and a confirming en¬ 
ergy deposit in an electromagnetic calorimeter (EMCal). 
The RICH uses CO 2 gas at atmospheric pressure as a 
Cerenkov radiator. Electrons and pions begin to ra¬ 
diate in the RICH at px > 20 MeV/c and px > 4.9 
GeV/c, respectively. The EMCal is composed of four 
sectors in each arm. The bottom two sectors of the east 







5 


arm are lead-glass and the other six are lead-scintillator. 
The energy resolution of the EMCal is ctb/E ~ 4.5% © 
8.3/v'£;(GeV) and ctb/E ~ 4.3% © l.l/y^E{GeV) for 
lead-scintillator and lead-glass, respectively. 



FIG. 1. (Color Online) (a) A schematic view of the PHENIX 
detector configuration for the 2011 run. (b) A schematic view 
of the VTX detector with the individual ladders shown. 


C. The VTX detector 

In 2011, the central detector was upgraded with the 
VTX detector as shown in Fig. In addition, a new 
beryllium beam pipe with 2.16 cm inner diameter and 760 
/rm nominal thickness was installed to reduce multiple¬ 
scattering before the VTX detector. 

The VTX detector [HHtiin] consists of four radial layers 
of silicon detectors as shown in Fig. [^b). The detector 
is separated into two arms, each with nominal accep¬ 
tance A(/) « O.Stt centered on the acceptance of the outer 
PHENIX central arm spectrometers. The detector cov¬ 
ers pseudorapidity \r]\ < 1.2 for collisions taking place at 
z = 0. The VTX can precisely measure the vertex posi¬ 
tion of a collision within \z\ < 10 cm range of the center 


of the VTX. 

The two inner layers, referred to as BO and Bl, of the 
VTX detector comprise silicon pixel detectors, as detailed 
in Ref. [41]. BO (Bl) comprises 10 (20) ladders with a 
central radial position of 2.6 (5.1) cm. The silicon pixel 
technology is based on the ALICElLHCb sensor-readout 
chip [42], which was developed at CERN. Each ladder 
is electrically divided into two independent half-ladders. 
Each ladder comprises four sensor modules mounted on 
a mechanical support made from carbon-fiber compos¬ 
ite. Each sensor module comprises a silicon pixel sensor 
with a pixel size of 50 /rm(())) x 425 /im(z) bump-bonded 
with four pixel readout chips. One pixel readout chip 
reads 256 (0)x 32 {z)= 8192 pixels and covers approxi¬ 
mately 1.3 cm 1.4 cm (Az) of the active area of 

the sensor. The position resolution is = 14.4 /im in 
the azimuthal direction. 

The two outer layers of the VTX detector, referred to 
as B2 and B3, are constructed using silicon stripixel sen¬ 
sors, as detailed in Ref. [41]. The B2 (B3) layer comprises 
16 (24) silicon stripixel ladders at a central radial distance 
of 11.8 (16.7) cm. The stripixel sensor is a novel silicon 
sensor, and is a single-sided, N-type, DC-coupled, two- 
dimensional (2-D) sensitive detector |43l|44|. One sensor 
has an active area of approximately 30 mm x 60 mm, 
which is divided into two independent sectors of 30 mm 
X 30 mm. Each sector is divided into 384 x 30 pixels. 
Each pixel has an effective size of 80 /rm {(j)) x 1000 /xm 
(z), leading to a position resolution of <j^—2'i ^m. A pixel 
comprises two implants (A and B) interleaved such that 
each of the implants registers half of the charge deposited 
by ionizing particles. There are 30 A implants along the 
beam direction, connected to form a 30 mm long X-strip, 
and 30 B implants are connected with a stereo angle of 
80 mrad to form a U-strip. X-strip and U-strip are visu¬ 
alized in [33]. When a charged particle hits a pixel, both 
the X- and the U-strip sharing the pixel register a hit. 
Thus the hit pixel is determined as the intersection of 
the two strips. The stripixel sensor is read out with the 
SVX4 chip developed by a FNAL-LBNL Collaboration 

03 ' 

The total number of channels in the VTX pixel and 
stripixel layers is 3.9 million pixels and 0.34 million strips. 
The compositions of the pixel and strip are illustrated in 
[311 mi- The main characteristics of the VTX detector 
are summarized in Table [H 


III. ANALYSIS 
A. Overview 

The purpose of the analysis is to separate the electrons 
from charm and bottom hadron decays. The life time of 
B mesons (ctbo= 455 /xm, CTg± = 491 /xm [36]) is sub¬ 
stantially longer than that of D mesons {ergo = 123 /xm, 
CTg± = 312 /xm) and the decay kinematics are differ¬ 
ent. This means that the distribution of values for the 
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TABLE I. A summary of the VTX detector. For each layer (BO to B3), the detector type, the central radius (r), ladder length 
(Z), sensor thickness (t), sensor active area {Acp x Az), the number of sensors per ladder {Ns), the number of ladders (Nl), 
pixel/strip size in <j) (Acj)) and z {Az), the nnmber of read-out channels (Nch), and the average radiation length including the 
support and on-board electronics (Ao) are given. 



type 

r(cm) 

Z(cm) 

t (pm) 

sensor active area 
A(()(cm) Az(cm) 

Ns 

Nl 

pixel/strip 
A(f) (pm) 

size 

Az (/^m) 

Ngh 

Vo(%) 

BO 

pixel 

2.6 

22.8 

200 

1.28 

5.56 

4 

10 

50 

425 

1.3 X 10® 

1.3 

B1 

pixel 

5.1 

22.8 

200 

1.28 

5.56 

4 

20 

50 

425 

2.6 X 10® 

1.3 

B2 

stripixel 

11.8 

31.8 

625 

3.07 

6.00 

5 

16 

80 

3 X lO"' 

1.2 X 10® 

5.2 

B3 

stripixel 

16.7 

38.2 

625 

3.07 

6.00 

6 

24 

80 

3 X 10"^ 

2.2 X 10® 

5.2 


distance of closest approach (DCA) of the track to the 
primary vertex for electrons from bottom decays will be 
broader than that of electrons from charm decays. There 
are other sources of electrons, namely Dalitz decays of 7r° 
and T], photon conversions, decays, and J/tp —)■ e+e” 
decays. With the exception of electrons from K^d, decays, 
these background components have DCA distributions 
narrower than those from charm decay electrons. Thus 
we can separate 6 —>■ e, c —?► e and background electrons 
via precise measurement of the DCA distribution. 


In the first step of the analysis, we select good events 
where the collision vertex is within the acceptance of the 
VTX detector, and its function is normal (Sec. IllB). We 
then reconstruct electrons in the PHENIX central arms 
(Sec. IIIC). The electron tracks are then associated with 
hits in the VTX detector and their DCA is measured 
(Sec. HID). At this point we have the DCA distribution 
of inclusive electrons that has contributions from heavy 
flavor (& —7> e and c e) and several background compo¬ 
nents. 


The next step is to determine the DCA shape and nor¬ 


malization of all background components (Sec. IIIEl. 
They include mis-identified hadrons, background elec¬ 
trons with large DCA caused by high-multiplicity effects, 
photonic electrons (Dalitz decay electrons, photon con¬ 
versions), and electrons from and quarkonia decays. 
The shapes of the DCA distributions of the various back¬ 
ground electrons are determined via data driven methods 
or Monte Carlo simulation. We then determine the nor¬ 
malization of those background electron components in 
the data (Sec. IIIF). 


Because the amount of the VTX detector material 
is substantial (13% of one radiation length) the largest 
source of background electrons is photon conversion 
within the VTX. We suppres s this background by a con¬ 
version veto cut (Sec. HIE3) 


Once the shape and the normalization of all back¬ 
ground components are determined and subtracted, we 
arrive at the DCA distribution of heavy flavor decay elec¬ 
trons that can be described as a sum of & —?► e and c —?► e 
DCA distributions. The heavy flavor DCA distribution 
is decomposed by an unfolding method (Sec. IIIC). 


B. Event selection 

The data set presented in this analysis is from Au-|-Au 
collisions at ,/s^ = 200 GeV recorded in 2011 after the 
successful commissioning of the VTX detector. As de¬ 
tailed earlier, the MB Au-|-Au data sample was recorded 
using the BBC trigger sampling 96 ± 3% of the inelastic 
Au-|-Au cross section. A number of offline cuts were ap¬ 
plied for optimizing the detector acceptance uniformity 
and data quality as described below. After all cuts, a 
data sample of 2.4x10® Au-I-Au events was analyzed. 


1. z-vertex selection 

The acceptance of the PHENIX central arm spectrom¬ 
eters covers collisions with z-vertex within ± 30 cm of the 
nominal interaction point. The VTX detector is more 
restricted in |z| acceptance, as the BO and B1 layers 
cover only \z\ < 11.4 cm. Thus the BBC trigger se¬ 
lected only events within the narrower vertex range of 
I-ZbbcI < 12 cm. In the offline reconstruction, the tracks 
reconstructed from VTX information alone are used to 
reconstruct the Au-|-Au collision vertex with resolution 
= 75 /im. All Au-f Au events in the analysis are re¬ 
quired to have a z-vertex within ±10 cm as reconstructed 
by the VTX. 


2. Data quality assurance 

Due to a number of detector commissioning issues in 
this first data taking period for the VTX, the data qual¬ 
ity varies substantially. Therefore we divide the entire 
2011 Au±Au data taking period into four periods. The 
acceptance of the detector changes significantly between 
these periods. 

In addition, several cuts are applied to ensure the qual¬ 
ity and the stability of the data. Applying electron iden¬ 
tification cuts described in Sec. EH C 2[ the electron to 
hadron ratios were checked for each run, a continuous 
data taking period typically lasting of order one hour, 
and three runs out of 547 with ratios outside of 5a from 
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the mean were discarded. The B2 and B3 stripixel lay¬ 
ers had an issue in stability of read-out electronics where 
some of the sensor modules would drop out, resulting 
in a reduced acceptance within a given run. Additional 
instabilities also existed in the BO and B1 pixel layers. 
Detailed channel by channel maps characterizing dead, 
hot, and unstable channels were generated for all layers 
within a given run. These maps were used to mask dead, 
hot, and unstable channels from the analysis, as well as 
to define the fiducial area of the VTX in simulations. 

During this first year of data taking, the instability 
of the read-out electronics discussed above caused sig¬ 
nificant run-to-run variations in the acceptance and ef¬ 
ficiency of the detector. It is therefore not possible to 
reliably calculate the absolute acceptance and efficiency 
correction while maintaining a large fraction of the to¬ 
tal data set statistics. Instead, we report on the relative 
yields of charm and bottom to total heavy flavor. We 
have checked that the DCA distributions are consistent 
between running periods and are not impacted by the 
changing acceptance. Thus we can measure the shape of 
the DCA distribution using the entire data set. In the 
following, we use the shape of the measured DCA distri¬ 
bution only to separate 6 —> e and c —>■ e components. 

C. Electron reconstruction in central arms 

1. Track reconstruction 

Charged particle tracks are reconstructed using the 
outer central arm detectors, DC and PC, as detailed in 
Ref. [IS] . The DC has six types of wire modules stacked 
radially, named XI, Ul, VI, X2, U2, and V2. The X 
wires run parallel to the beam axis in order to measure 
the 0-coordinate of the track and the U and V wires have 
stereo angles varying from 5.4 to 6.0 degrees. Tracks are 
required to have hits in both the XI and X2 sections 
along with uniquely associated hits in the U or V stereo 
wires and at least one matching PC hit, to reduce mis- 
reconstructed tracks. The track momentum vector is de¬ 
termined assuming the particle originated at the Au-|-Au 
collision vertex as reconstructed by the BBC. 

2. Electron identification 

Electron candidates are selected by matching tracks 
with hits in the RICH and energy clusters in the EMCal. 
The details on the electron selection cuts are given in 
Ref. [H]. In this analysis we select electron candidates 
within 1.5 < pt [GeV/c] < 5.0, and we briefly describe 
the cuts in the RICH and EMCal below. 

Cerenkov photons from an electron track produce a 
ring-shaped cluster in the RICH. At least three associ¬ 
ated PMT hits are required in the RICH and a ring-shape 
cut is applied. The center of the ring is required to be 
within 5 cm of the track projection. The probability 


that the associated cluster in the EMCal comes from an 
electromagnetic shower is calculated based on the shower 
shape. Based on that probability, tracks are selected in a 
way that maintains high efficiency for electrons while re¬ 
jecting hadrons. Further, the energy (E) in the EMCal is 
required to match the track determined momentum (p). 
This match is calculated as dep = {E/p — PEjr^l'^Elrn 
where and cr^/p are the mean and standard devia¬ 
tion respectively of a Gaussian fit to the E jp distribution, 
determined as a function of momentum (see Fig. [^. A 
cut of dep > — 2 is used to further reject hadrons that 
have dJi E/p ratio < I, because they do not deposit their 
full energy in the EMCal. 

In high-multiplicity Au-l-Au events there is a signifi¬ 
cant probability for a random association between the 
track and hits in the RICH and EMCal. This mis- 
identified hadron probability is estimated as follows. The 
z < 0 and z > 0 sides of the RICH have their hits 
swapped in software, and the tracks are re-associated 
with RICH hits. Because the two longitudinal sides of 
the RICH are identical, this gives a good estimate of the 
random hadron background in the electron sample. 

The distribution of electron candidates at px =2.0-2.5 
GeV/c for the normalized EMCal energy to track mo¬ 
mentum ratio, dep defined above, is shown in Fig. 
There is a large peak near zero from true electrons as 
expected and a clear low-side tail from mis-identified 
hadron. Also shown is the result of the above swap 
method. The difference between the data and the “swap” 
distribution (red) is explained as contributions from off- 
vertex electrons caused by conversions from the outer 
layer of the VTX and weak decay. In the final account¬ 
ing for all contributions to the identified-electron DCA 
distribution, we utilize this swap method to statistically 
estimate the contribution of mis-identified hadron in each 
Pt selection as detailed in Section [IH E 1[ 


2<p^<2.5 



FIG. 2. (Color Online) Matching variable between the re¬ 
constructed track momentum (p) and the energy measured in 
the EMCal (A): dep = {E/p — PE/p)/<^E/r>- The black dis¬ 
tribution is for identified electrons with pr = 2.0-2.5 GeV/c, 
and the red distribution is the estimated contribution from 
mis-identified electrons via the RICH swap-method. 






D. DCA measurement with the VTX 


3. The primary vertex reconstruction 


Charged particle tracks reconstructed in the central 
arms must be associated with VTX hits in order to cal¬ 
culate their DCA. Three-dimensional (3-D) hit positions 
in the 4 layers of VTX are reconstructed. For each col¬ 
lision, the primary vertex is reconstructed by the VTX. 
Then central arm tracks are associated with hits in the 
VTX, and VTX-associated tracks are formed. Finally, 
the DCA between the primary vertex and the VTX- 
associated tracks are measured. 


1. VTX alignment 

In order to achieve good DCA resolution to separate 
6 —?► e and c —>■ e, alignment of the detector ladders 
to high precision is required. The detector alignment 
is accomplished via an iterative procedure of matching 
outer central arm tracks from the DC and PC to the 
VTX hits. The procedure is convergent for the position 
of each ladder. The alignment was repeated each time 
the detector was repositioned following a service access. 
The final alignment contribution to the DCA resolution 
in both (j) and z is a few tens of microns. 

8. VTX hit reconstruction 

For layers BO and Bl, clusters of hit pixels are formed 
by connecting contiguous hit pixels by a recursive cluster¬ 
ing algorithm. An average cluster size is 2.6 (6.7) pixels 
for the pixel (stripixel). The center of the cluster in the 
local 2-D coordinate system of the sensor is calculated as 
the hit position. 

For B2 and B3 layers, 2D hit points on the sensor are 
reconstructed from the X-view and the U-view. Hit lines 
in the X-view (U-view) are formed by clustering contigu¬ 
ous hit X-strips (U-strips) weighted by deposited charges, 
and then 2D hit points are formed as the intersections of 
all hit lines in X- and U- views. When one hit line in U- 
view crosses more than two hit lines in X-view, ghost hits 
can be formed, because which crossing point is the true 
hit is ambiguous. These ghost hits increase the number 
of reconstructed 2D hits approximately by 50% (30%) in 
B2 (B3) in central Au-|-Au collisions. The ghost hit rate 
was studied using a full GEANt3 simulation with the 
HIJING [IS] generator as input. However, because the 
occupancy of the detector at the reconstructed 2D hit 
point level is low, less than 0.1%, these ghost hits do not 
cause any signihcant issue in the analysis. 

The positions of all 2-D hits in the VTX are then trans¬ 
ferred into the global PHENIX 3-D coordinate system. 
Correction of the sensor position and orientation, deter¬ 
mined by the alignment procedure described in the previ¬ 
ous section, is applied in the coordinate transformation. 
The resulting 3-D hit positions in the global coordinate 
system are then used in the subsequent analysis. 


With the VTX hit information alone, charged particle 
tracks can be reconstructed only with modest momen¬ 
tum resolution 5p/p « 10% due to the limited magnetic 
field integrated over the VTX volume and the multiple 
scattering within the VTX. These tracks can be utilized 
to determine the collision vertex in three-dimensions (zg 
along the beam axis, and Xo-,yo in the transverse plane) 
for each Au-|-Au event under the safe assumption that 
the majority of particles originate at the collision vertex. 
This vertex position is called the primary vertex position. 

The position resolution of the primary vertex for each 
direction depends on the sensor pixel and strip sizes, the 
precision of the detector alignment, and the number of 
particles used for the primary vertex calculation and their 
momentum in each event. For MB Au-|-Au collisions, 
the resolution values are = 96 ^m, ay = 43 ^m, and 
az = 75 pm. The worse resolution in x compared to 
y is due to the orientation of the two VTX arms. For 
comparison, the beam profile in the transverse plane is 
^^umi « 90 pm in the 2011 Au-|-Au run. 

4- Association of a central arm track with VTX 

Each central arm track is projected from the DC 
through the magnetic field to the VTX detector. Hits 
in VTX are then associated with the track using a recur¬ 
sive windowing algorithm as follows. 

The association starts from layer B3. VTX hits in that 
layer that are within a certain (A^ x Az) window around 
the track projection are searched. If hits are found in this 
window, the track is connected to each of the found hits, 
and then projected inward to the next layer. In this case 
the search window in the next layer is decreased, because 
there is much less uncertainty in projection to the next 
layer. If no hit is found, the layer is skipped, and the 
track is projected inward to the next layer, keeping the 
size of the projection window. This process continues 
until the track reaches layer BO, and a chain of VTX hits 
that can be associated with the track is formed. The 
window sizes are momentum dependent and determined 
from a full GEANt3 simulation of the detector so that the 
inefficiency of track reconstruction due to the window size 
is negligible. 

After all possible chains of VTX hits that can be as¬ 
sociated with a given central arm track are found by the 
recursive algorithm, a track model ht is performed for 
each of these possible chains, and the of the fit, 
is calculated. The effect of multiple scattering in each 
VTX layer is taken into account in calculation of Xvtx- 
Then the best chain is chosen based on the value of Xvtx 
and the number of associated hits. This best chain and 
its track model are called a VTX-associated track. Note 
that at most one VTX-associated track is formed from 
each central arm track. 

In this analysis we require that VTX-associated tracks 
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have associated hits in at least the first three layers, i.e. 
BO, Bl, and B2. An additional track requirement is 
X?tx/NDF < 2 for pt < 2 GeV/c and Xvtx/NDF < 3 
for pt > 2 GeV/c, where NDF is the number of degrees 
of freedom in the track fit. 


5. DCAt and DCAi 


Using the primary vertex position determined above, 
the DGA of a track is calculated separately in the trans¬ 
verse plane (DCA-r) and along the beam axis (DGA^). 
Because by design the DCAt has a better resolution than 
DCAt, we first find DCAt with a track model of a circle 
trajectory assuming the uniform magnetic field over the 
VTX. We define DCAt as 


DCAt = L-R, (1) 


where L is the distance from the collision vertex to the 
center of the circle defining the particle trajectory, and 
R is the radius of the circle as shown in Fig. DCAt is 
the distance between the z-coordinate of the point DCAt 
found and z-coordinate of the primary vertex. 


It is notable that DCAt has a sign in this defini¬ 
tion. The distinction between positive and negative val¬ 
ues of DCAt— whether the trajectory is bending towards 
or away from the primary vertex—is useful since cer¬ 
tain background contributions have asymmetric distri¬ 
butions in positive and negative DCAt, as discussed in 
section HIE For electrons, the positive side of DCAt 
distribution has less background contribution. There is 
no such positive/negative asymmetry in DCAt. 





FIG. 3. (Color Online) Illustration of the definition of DCAt 
= L - R in the transverse plane. 


6. DCA measurement 


For each VTX-associated track, the DCA is calcu¬ 
lated separately in the radial and longitudinal direction 
(DCAt and DCAt) from the track model and the pri¬ 
mary vertex position. Shown in Fig. is the resulting 
DCAt and DCAt distributions for all VTX-associated 
tracks withpT = 2.0-2.5 GeV/c. Since the vast majority 
of charged tracks are hadrons originating at the primary 
vertex, we observe a large peak around DCAt, DCAt = 
0 that is well fit to a Gaussian distribution where the a 
represents the DCAt, DCAt resolution. A selection of 
I DCAt I < 0.1 cm is applied to reduce background. 

There are broad tails for |DCAt | > 0.03 cm. Monte 
Carlo simulation shows that the main source of the broad 
tails is the decay of long lived light hadrons such as A and 


Kl 

The DCAt resolution as a function of the track pT is 
extracted using a Gaussian fit to the peak and is shown in 
Fig. 0 c). The DCAt resolution is approximately 75 pm 
for the 1.0-1.5 GeV/c bin and decreases with increasing 
Pt as the effect of multiple scattering becomes smaller for 
higher pr- The DCAt resolution becomes less than 60 
pm for Pt > 4 GeV/c, where it is limited by the position 
resolution of the primary vertex. 

We divide the electrons into five pr bins and show the 
DCAt distributions for each in Fig. These distribu¬ 
tions are in integer-value counts and are not corrected 
for acceptance and efficiency. The DCA distributions in¬ 
clude various background components other than heavy 
flavor contributions. The background components are 
also shown in the figure and are discussed in the next 
section (Section IIIEl. 

While the DCAt distributions in Fig. are plotted 
within |DCAt| < 0.15 cm, only a |DCAt| < 0.1 cm is 
used in the analysis to extract the charm and bottom 
yield described later. At large DCAt, the distribution is 
dominated by high-multiplicity background (Sec. HIE2) 
and therefore provides little constraint in the extraction 
of the charm and bottom contributions. 


E. DCA distribution of Background Components 

The sample of candidate electron tracks that pass all 
the analysis cuts described above contains contributions 
from a number of sources other than the desired elec¬ 
trons from semi-leptonic decays of charm and bottom 
hadrons. In order to extract the heavy flavor contri¬ 
butions, all background components must be fully ac¬ 
counted for and their DCAt shapes as a function of pt 
incorporated. These background components are listed 
in the order presented below. 

1. Misidentified hadrons 

2. High-multiplicity background 

3. Photonic electrons 
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FIG. 4. Distance-of-closest-approach distributions for (a) 
along the beam axis DCA_l and (b) transverse plane DC At for 
all VTX-associated tracks in Au+Au at = 200 GeV in 

the range 2.0 < pr [ GeV/c] < 2.5. (c) The DCAt resolution 
as a function of pr for all tracks. 


4. Kaon decay electrons 

5. Heavy-quarkonia decay electrons 

As described in this and the following section, all back¬ 
ground components are constrained by PHENIX mea¬ 
surements in Au-bAu and are fully simulated through a 
GEANt3 description of the detector. This method is sim¬ 
ilar to the cocktail method of background subtraction 


used in the previous analysis of inclusive heavy flavor 
electrons [l2]. 

Next, we describe these background sources and their 
DCA distributions. The first two components are caused 
by detector and multiplicity effects. DCA distributions 
and normalization of these two components are deter¬ 
mined by data driven methods, as detailed in this sec¬ 
tion. The last three components are background elec¬ 
trons that are not the result of semi-leptonic decays of 
heavy flavor hadrons. Their DCA distributions are de¬ 
termined by Monte Carlo simulation, and their normal¬ 
ization is determined by a bootstrap method described 
in section III F Of those background electrons, photonic 
electrons are the dominant contribution. We dev eloped a 
conversion veto cut to suppress this background (III E 31. 


1. Mis-identified hadron 


As detailed in the discussion on electron identification, 
there is a nonzero contribution from mis-identified elec¬ 
trons. This contribution is modeled via the RICH swap- 
method described in Section El C 2[ From this swap 
method, we obtain the probability that a charged hadron 
is mis-identified as an electron as a function of pp ■ This 
probability is then applied to the DCA distribution of 
charged hadrons to obtain the DCA distribution of mis- 
identified hadrons. 

The resulting DCAt distribution is shown in each 
panel of Fig. Note that this component is properly 
normalized automatically. For each pp bin, the DCA 
distribution of mis-identified prompt hadrons has a nar¬ 
row Gaussian peak at DCAt = 0. The broad tails for 
large |DCAt | are mainly caused by decays of A and K^. 
In all pp bins the magnitude of this background is no 
more than 10% of the data for all DCAt 


2. High-multiplicity background 

Due to the high multiplicity in Au-|-Au collisions, an 
electron candidate track in the central arms can be asso¬ 
ciated with random VTX hits. Such random associations 
can cause a background that has a very broad DCAt dis¬ 
tribution. Although the total yield of this background is 
only ~ 0.1% of the data, its contribution is significant at 
large DCAt where we separate 5 —?> e and c —>■ e. 

To evaluate the effect of event multiplicity on the 
reconstruction performance, we embed simulated sin¬ 
gle electrons— i.e. the response of the PHENIX detec¬ 
tor to single electrons that is obtained from a GEANt3 
simulation—into data events containing VTX detector 
hits from real Au-I-Au collisions. The events are then pro¬ 
cessed through the standard reconstruction software to 
evaluate the reconstruction performance in MB Au-|-Au 
collisions. 

The reconstructed DCAt and DCA^ for embedded 
primary electrons in MB Au-l-Au collisions is shown in 
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FIG. 5. (Color Online) DCAt distributions for electrons in MB Au+Au at = 200 GeV that pass the reconstruction and 

conversion veto cut in the indicated five electron-pT selections. Also shown are the normalized contributions for the various 
background components detailed in Section [ill E| 


Fig. § Here the histograms, labeled as “Single Elec- trons”, show the reconstructed DCAp and DCAp dis- 
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DCA-p [cm]" 



DCA^I [cm] 


FIG. 6. (Color Online) Simulated primary electron (a) 
DCAt and (b) DCA_l distribution before and after embed¬ 
ding in real Au-|-Au data. 


butions have broad tails at large |DCAr| and |DCAi|. 

As shown in Fig. [^b), tracks with |DCAl| > 0.13 cm 
are dominated by random associations, as they are not 
present in the “Single Electron” sample. We therefore 
use the DCAt distribution for tracks with large |DCAt| 
as an estimate of this random high-multiplicity back¬ 
ground. We choose the region 0.13 < |DCAt| cm < 0.18 
to represent this background, and restrict our signal to 
|DCAt| <0.1 cm. The DCA-p distribution of tracks 
with 0.13 < |DCAt| cm < 0.18 must be normalized in 
order to be used as an estimate of the high-multiplicity 
background for tracks within |DCAt| <0.1 cm. This 
normalization is determined by matching the integrated 
yield of embedded primary electrons in each |DCAt| re¬ 
gion for 0.08 < DCAt cm < 0.2, as shown in the inlay 
of Fig. [^b). The region 0.08 < DCAt cm < 0.2 is dom¬ 
inated by random associations, as shown in Fig. §a), 
and is therefore safe to use for determining the normal¬ 
ization. The normalization of the high-multiplicity back¬ 
ground is determined to be 2.89 ± 0.29. The red filled 
circles in Fig. I^a) show the embedded DCAt distri¬ 
bution with large DCAt (0.13 < |DCAt| cm < 0.18). 
This distribution agrees with the embedded DCAt dis¬ 
tribution (blue filled triangles in Fig. for large DCAt. 
This demonstrates that the tails for large DCAt are well 
normalized by the distribution of electrons with large 
DCAt. However, there is a small excess in the region 
0.05 < |DCAt| cm < 0.10 that is not accounted for by 
the distribution with large DCAt. We address this excess 
in the systematic uncertainties, as described in Sec. [Ill H] 
where it is found to have only a small effect on the ex¬ 
traction of 6 —>■ e and c —>■ e. 

In each panel of Fig. [^ the high-multiplicity back¬ 
ground is shown as a red line. It is determined from 
the DCAt distribution of the data within 0.13 < 
|DCAt| cm < 0.18, as described above. The number 
of electron tracks in the large DCAt region is small. We 
therefore fit the resulting DCAt data in each pT bin with 
a smooth function to obtain the shape of the red curves 
shown in Fig. [^ A second order polynomial is used in 
the lowest pt bin, where there are enough statistics to 
constrain it. The higher px bins are fit with a constant 
value. All curves are multiplied by the same normaliza¬ 
tion factor, determined from embedded simulations as 
described above. 


tributions of primary electrons before embedding. The 
DCAt distribution comprises a narrow Gaussian with 
no large DCAt tail and the DCAt distribution com¬ 
prises a similar, but slightly broader, Gaussian with no 
large tail. The blue filled triangles show the DCAt 
and DCAt distributions after embedding. The DCAt 
and DCAt distributions comprise a Gaussian peaked at 
DCAt(DCAt) 0 which is consistent with the distri¬ 
bution before embedding. This demonstrates that the 
DCA resolution of the VTX is not affected by the high 
multiplicity environment. However, the embedded distri- 


3. Photonic electrons and conversion veto cut 

Photon conversions and Dalitz decays of light neutral 
mesons (tt*^ and rj) are the largest electron background. 
We refer to this background as photonic electron back¬ 
ground as it is produced by external or internal conver¬ 
sion of photons. 

The PHENIX Collaboration has previously published 
the yields of 7r° and rj mesons in Au-|-Au collisions at 
= 200 GeV [1^ [SD]. In addition to the electrons 
from Dalitz decays of these mesons, the decay photons 
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may convert to an e“'"e“ pair in the detector material in 
the beam pipe or each layer of the VTX. The PHENIX 
Collaboration has also published the yields of direct pho¬ 
tons in Au-I-Au collisions at y/s~^ = 200 GeV [31 I3T] . 
that can also be a source for conversions. 

In principle with these measured yields, combined with 
simple decay kinematics and a detailed GEANt3 descrip¬ 
tion of the detector material and reconstruction algo¬ 
rithm, one could fully account for these photonic electron 
contributions as a function of DCA-r and pt- However, 
systematic uncertainties on the measured yields for the 
TT^, ry, and direct photons would then dominate the uncer¬ 
tainty of the heavy flavor electron extraction. Therefore, 
we utilize the VTX detector itself to help reject these 
contributions in a controlled manner. 

We require that at least the first three layers of the 
VTX have hits associated with the electron track. Con¬ 
versions in B1 and subsequent layers are rejected by the 
requirement of a BO hit, leaving only conversions in BO 
and the beam pipe. The requirement of B1 and B2 hits 
enables us to impose a conversion veto cut, described be¬ 
low, that suppresses conversions from the beam pipe and 
BO. 

The conversion veto cut rejects tracks with another 
VTX hit within a certain window in A</> and Az around 
hits associated with a VTX-associated track. Photons 
that convert to an e“'“e“ pair in the beam pipe will leave 
two nearby hits in the first layer (BO) and/or subsequent 
layers of the VTX, and thus be rejected by the conversion 
veto cut. Similarly, conversions in BO will result in two 
nearby hits in the second layer (Bl) and/or subsequent 
outer layers. The same is true for e+e“ from a Dalitz 
decay, though with a larger separation due to a larger 
opening angle of the pair. 

Figurej^a) shows distribution of chrg Afp of hits in BO 
relative to the electron track, where chrg is the charge of 
the track. The red (circle) histogram shows the data in 
MB Au-|-Au collisions. If the track at the origin is not an 
electron, we have a flat distribution due to random hits 
in the detector. These random hits have been subtracted 
in Fig.j^a). The transverse momentum of the electron 
track is in the interval 1 < pr GeV/c < 2. 

As mentioned above, these correlated hits around elec¬ 
tron tracks are caused by the partner e"*" or e~ of Dalitz 
decays or photon conversions. The left-right asymmetry 
of the distribution is caused by the fact that the part¬ 
ner track is separated from the electron track by the 
magnetic field and the direction of the separation is de¬ 
termined by the charge of the electron track. In the dis¬ 
tribution of chrg Ap, the partner track is bent towards 
the positive direction. 

The black (triangle) histogram in Fig. [^a) shows the 
distribution from Monte Carlo simulations. In the simu¬ 
lation, the response of the PHENIX detector to single 7r°s 
is modeled by GEANt3, and the resulting hits in the VTX 
and the central arms are then reconstructed by the same 
reconstruction code as the data. The correlated hits in 
the simulation are caused by the Dalitz decay of and 


chrg A(|) at BO 



Veto residual BO ( chrg A(|) vs ) 



FIG. 7. (Color Online) (a) Distribution of correlated hits 
in BO near electron tracks for 1 < pt < 2 GeV/c. The red 
(circle) points are from Au+Au data and the black (triangle) 
points are from Monte Carlo simulation. The insert in (a) 
illustrates the electron pairs from Dalitz decays, (b) The win¬ 
dow of the conversion veto cut for BO layer (hatched) and the 
hit distribution near electron track in 2D space of chrg A(j> vs 
Pt of electrons in Au-|-Au collisions. (See the text for details). 


photon conversion in the material of the beam pipe and 
the VTX itself. The simulation reproduces the data well 
for chrg A(j> > 0. There is a difference between the data 
and the simulation for chrg Ap < 0. This is caused by a 
subtle interplay between the conversions and high mul¬ 
tiplicity effects. The difference disappears for peripheral 
collisions. Similar correlated hits are observed in BI to 
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B3 layers in the data and they are also well explained by 
the simulation. 

We define a “window” of the conversion veto cut 
around an electron track in each layer BO to B3 and 
require that there is no hit other than the hit associ¬ 
ated with the electron track in the window. Since a pho¬ 
tonic electron (Dalitz and conversion) tends to have a 
correlated hit in the window, as one can see in Fig. 
this conversion veto cut rejects photonic background. A 
larger window size can reject photonic background more 
effectively, but this can also reduce the efficiency for the 
heavy flavor electron signal due to random hits in the 
window. The window for the conversion veto cut is a 
compromise in terms of the rejection factor on photonic 
backgrounds and efficiency for heavy flavor electrons. We 
optimized the size of the window of the conversion veto 
cut based on a full GEANt3 simulation. 

The red hatched area shown in Fig. shows the win¬ 
dow of the conversion veto cut in layer BO. The window 
size is asymmetric since correlated hits are mainly in the 
positive side of chrg Acj). The window size is reduced for 
higher electron px since the distribution of correlated hits 
becomes narrower for higher px- The windows for B1-B3 
are similarly determined based on GEANt3 simulation. 

Figure shows the survival fraction of the conversion 
veto cut for electrons from photon conversions and Dalitz 
decays as a function of electron px from a full GEANt3 
simulation of the detector with hits run through the re¬ 
construction software. The survival probability for con¬ 
versions is less than 30% at pr = 1 GeV/c and decreases 
further at higher px- The survival probability for Dalitz 
decays is higher since a Dalitz decay partner is more likely 
to fall outside of the window of the conversion veto cut 
due to the larger opening angle. Also shown in Fig. [^is 
the survival fraction of electrons from heavy flavor decays 
which pass the conversion veto cut (5'hf)- As expected, 
their efficiency for passing the conversion veto cut is quite 
high and px independent. 

The efficiencies shown in Fig. are calculated without 
the Au+Au high-multiplicity that may randomly pro¬ 
vide a hit satisfying the conversion veto cut. Since these 
are random coincidences, they are a common reduction 
for all sources including the desired signal — heavy fla¬ 
vor electrons. This common reduction factor, 5randorm is 
measured from the reduction of the hadron track yield 
by the conversion veto cut to be ~ 35% at pt = 1 GeV/c 
to ~ 25% at pt = 5 GeV/c for MB Au+Au collisions. 
Note that when we determine the DCAt distribution of 
the various background components using a full GEANt3 
simulation we apply the same conversion veto cuts. 

The DGAt distributions from photonic background 
processes that survive the conversion veto cut are shown 
in Fig. The means of the DGA^ distributions from 
Dalitz decays and conversions are shifted to negative 
DGAt values due to the mis-reconstruction of the mo¬ 
mentum caused by the assumption that the tracks orig¬ 
inate at the primary vertex, as explained in the next 
paragraph. The shift is largest at the lowest px bin and 



FIG. 8. (Color Online) The survival rate as a function of 
electron px (px) for electrons from photon conversion (black), 
Dalitz decay of (red), rj (green), electrons from direct pho¬ 
ton (blue) and heavy flavor decay electrons (dark orange). 


decreases with increasing px ■ 

For Dalitz electrons, the shift is due to the energy loss 
via induced radiation (bremsstrahlung). The total radi¬ 
ation length of the VTX is approximately 13% as shown 
in Table Thus a Dalitz electron coming from the pri¬ 
mary vertex loses approximately 1 — « 12% of its 

energy on average when it passes through the VTX. The 
momentum measured by the DC is close to the one af¬ 
ter the energy loss due to the reconstruction algorithm. 
Since the momentum determined by the DC is used when 
projecting inward from the hit in BO to the primary ver¬ 
tex and in calculation of DCAt, this results in a slight 
shift in the DCAt distribution. This effect is fully ac¬ 
counted for in the DCAt template of Dalitz electrons 
since it is generated through the full GEANt3 and recon¬ 
struction simulation. 

In the case of conversions, the effect is even larger, 
as one can clearly see in Fig. While a photon goes 
straight from the primary vertex to the beam pipe or 
BO layer where it converts, DCAt is calculated assum¬ 
ing that the electron track is bent by the magnetic field. 
Thus the DCAt distribution is shifted by the difference 
of the actual straight line trajectory and the calculated 
bent trajectory. Again, this is fully accounted for with 
the full GEANt3 simulation. The effect is verified by se¬ 
lecting conversion electrons with a reversed conversion 
veto cut. 
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A- 

The background from i^es decays {Kg, ^ ei^n) 

contributes electrons over a broad range of DC At due 
to the long lifetime of the kaons. Both contributions are 
determined using pythia and a full GEANtS simulation, 
taking into account the exact track reconstruction, elec¬ 
tron identification cuts, and conversion veto cut. The re¬ 
sulting DCAt distribution for these kaon decays is shown 
in Fig. As expected, though the overall yield is small, 
this contributes at large DCAt in the lower pr bins and 
is negligible at higher px- 


5. Quarkonia 

Quarkonia (T/V' and T) decay into electron pairs. Due 
to the short lifetime, these decays contribute to electrons 
emanating from the primary vertex. The J/^/> yields in 
Au-I-Au collisions at = 200 GeV have been mea¬ 

sured by the PHENIX Collaboration [^. The detailed 
modeling of these contributions out to high px is detailed 
in Ref. [T^. While these measurements include a small 
fraction of R —>■ J/ijj decays, all J/^’s are considered 
prompt when modeling the DCAt distribution. The J/ip 
contribution is shown in Fig. and is quite small and 
peaked about DCAt = 0 as expected. Thus, the system¬ 
atic uncertainty from the quarkonium yields in Au-|-Au 
collisions is negligible in all electron px bins. 


Some backgrounds have the same ratio to signal regard¬ 
less of the year the data was collected, while others will 
differ due to the additional detector material added by 
the VTX. 

The invariant yield in Au-f Au collisions at 
= 200 GeV of heavy flavor electrons and back¬ 
ground electrons from Dalitz decays is a physical 
observable independent of the year the data was taken. 
Thus we can use the ratio of heavy flavor/Dalitz that is 
determined in the 2004 analysis in the 2011 data. On the 
other hand, the invariant yield of conversion electrons 
depends on the detector material present and is thus 
different in the 2011 data taking period with the VTX 
installed compared with the 2004 data. We account for 
this difference by calculating the fraction of nonphotonic 
electrons in the 2011 data. A detailed description of the 
normalization procedure is given in Appendix |VI| 

With this bootstrapped normalization completed, the 
correctly normalized background components are shown 
for all five px bins vs DGAt in Fig. Note that the 
normalization of mis-identified hadron and random back¬ 
ground is determined from the data as explained in sec¬ 
tions |1II E 1 and IIIE2[ respectively. The electron yield 
beyond the sum of these background components is from 
the combination of charm and bottom heavy flavor elec¬ 
trons. 


G. Unfolding 


F. Normalization of electron background 
components 

If the detector performance were stable, we could con¬ 
vert the DCAt distributions from counts into absolutely 
normalized yields. Then one could straightforwardly 
subtract the similarly absolutely normalized background 
contributions described above—with the normalization 
constrained by the previously published PHENIX yields 
for 7r°, r], etc. However, due to detector instability dur¬ 
ing the 2011 run, such absolute normalization of back¬ 
ground contributions can have a large systematic uncer¬ 
tainty. Thus we bootstrap the relative normalization of 
these background contributions utilizing our published 
Au-|-Au results [12] from data taken in 2004. 

The idea of the method is the following. PHENIX 
measured the invariant yield of open heavy flavor decay 
electrons from the 2004 dataset. In this 2004 analysis 
we first measured inclusive electrons (be. the sum of 
background electrons and heavy flavor electrons). We 
then determined and subtracted the background electron 
components from the inclusive electron yields to obtain 
the heavy flavor contribution. Thus the ratio of the back¬ 
ground components to the heavy flavor contribution were 
determined and published in [T2| . We use these ratios to 
determine the normalization of background components 
in the 2011 data, as described in the next paragraph. 


1. Introduction 

With the DCAt distributions as a function of electron 
Px and the various background components in hand, we 
proceed to extract the remaining charm and bottom com¬ 
ponents. If one knew the shape of the parent charm and 
bottom hadron px and rapidity distributions, one could 
calculate in advance the DCAt shape for electrons from 
each heavy flavor via a model of the decay kinematics. 
Since the decay lengths of charm and bottom hadrons are 
significantly different, they will yield different DCAt dis¬ 
tributions. In this case, one could simultaneously fit the 
DCAt distribution for each px bin with all background 
components fixed across px bins, and extract the one free 
parameter: the ratio of charm to bottom contributions. 
However, the px distribution of charm hadrons is known 
to be significantly modified in Au-|-Au collisions — see 
for example Ref. M- For bottom hadrons this is also 
likely to be the case. Therefore one does not know a pri¬ 
ori the heavy flavor DCAt distribution since it depends 
on the parent px distribution. 

Since the DCAt distributions for all electron px re¬ 
sult from the same parent charm and bottom hadron px 
spectrum, one can perform a simultaneous fit to all the 
electron px and DCAt data in order to hnd the most 
likely heavy flavor parent hadron px distributions. The 
estimation of a set of most likely model parameters using 
a simultaneous fit to data is often referred to as unfold- 
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ing. Statistical inference techniques are often employed 
to solve such problems; see for example the extraction of 
reconstructed jet cross sections [S3]. 

The DCA-r distributions are in counts and have not 
been corrected for the pT-dependent reconstruction effi¬ 
ciency in Au-|-Au collisions, and therefore hold no yield 
information. To further constrain the extraction of the 
charm and bottom components, we include the total 
heavy flavor electron invariant yield as measured by 
PHENIX [T2| in Au-|-Au collisions at ^/s^ = 200 GeV. 
This measurement is more accurate than currently avail¬ 
able with the 2011 data set, where the VTX acceptance 
changes with time. 

The unfolding procedure, us ing a p articular sampling 
method (described in Section III G 21, chooses a set of 
trial charm and bottom parent hadron yields. The trial 
set of yiel ds is mu ltiplied by a decay matrix (described 
in Section IIIG4), which encodes the probability for a 
hadron in a given px interval to decay to an electron at 
midrapidity as a function of electron px and DGA^. The 
resulting distributions of electron px and DCAy are com¬ 
pared with the measured data using a likelihood func¬ 


tion (described in Section III G 3). In order to dampen 
discontinuities and oscillatory beha vior, a p enalty upon 
the likelihood (described in Section [III G 5 ) is added to 
enforce smoothness in the resulting hadron px distribu¬ 
tions. 


denominator P(x) serves as an overall normalization of 
the combined likelihood P(x|0)7r(0) such that p(0|x) can 
be interpreted as a probability density. In this analysis, 
p(6|x) gives the probability for a set of charm and bottom 
hadron yields. 


e^{9,-,9b), (3) 

given the values of the measured electron data points 
X. Since we are only interested in the parameters which 
maximize p(0|x), we can dispense with the calculation of 
P(x), as it serves only as an overall normalization. 

Here 9 comprises 17 bins of both charm and bottom 
hadron px, yielding a 34-dimensional space which must 
be sampled from in order to evaluate p(0|x). To ac¬ 
complish this we employ a Markov Chain Monte Carlo 
(MCMC) algorithm to draw samples of 9 in proportion 
to p(0|x). This makes accurate sampling of multidimen¬ 
sional distributions far more efficient than uniform sam¬ 
pling. In implementation, it is in fact the right hand 
side of Eq. [^ that is sampled. The MCMC variant used 
here is an affine-invariant ensemble sampler described in 
Ref. [SS] and implemented as described in Ref. [33]. It 
is well suited to distributions that are highly anisotropic 
such as spectra which often vary over many orders of 
magnitude. 


2. Unfolding method 


3. Modeling the likelihood function 


Here we apply Bayesian inference techniques to the 
unfolding problem. A detailed pedagogical introduction 
to these techniques is given in Ref. [54] . Techniques in¬ 
volving maximum likelihood estimation or maximum a 
posteriori estimation, often used in frequentist statistics, 
can at best compute only a point estimate and confidence 
interval associated with individual model parameters. In 
contrast, Bayesian unfolding techniques have the impor¬ 
tant advantage of providing a joint probability density 
over the full set of model parameters. In this analysis, 
the vector of model parameters, 0, is the vector of parent 
charm and bottom hadron yields binned in px- 

Given a vector of measured data, x, and our vector of 
model parameters, 9, we use Bayes’ theorem 


p(0|x) = - 


P(x|0)^(0) 


P(x) 


( 2 ) 


to compute the posterior probability density p(0|x) 
from the likelihood P{x\9) and prior information 7r(0). 
The function P(x|0), quantifies the likelihood of observ¬ 
ing the data given a vector of model parameters. In 
frequentist statistics, the P(x|0) is often used alone to 
determine the best set of model parameters. Bayesian 
inference, on the other hand, allows for the inclusion of 
the analyzer’s a priori knowledge about the model pa¬ 
rameters, as encoded in 7r(0). The implementation of 
Tr{9) used in this analysis is discussed in Sec. HI G 5 The 


This analysis is based on 21 data points of total heavy 
flavor electron invariant yield, in the range 1.0- 

9.0 GeV/c from the 2004 data set [T3], and five electron 
DCAt distributions where j indexes each electron 

Px interval within the range 1.5-5.0 GeV/c from the 2011 
data set. Therefore, 

^ j-^data j^data j^data j^data j^data^ 

in Eq. [^ 

Our ultimate goal is to accurately approximate the 
posterior distribution over the parent hadron invariant 
yields 9 by sampling from it. Eor each trial set of hadron 
yields, the prediction in electron px, Y(0), and DCAt, 
D,(0), is calculated by 

Y{9) = MWec + MW6/j, (5) 

B,{9)=MP9, + Mp9b, ( 6 ) 

where and are decay matrices discussed in 

Section [ill G 41 We then evaluate the likelihood between 
the prediction and each measurement in the data sets 
Ydata |]^data|4^^^ jg customai'y, the logarithm of 
the likelihood function is used in practice. The combined 
(log) likelihood for the data is explicitly 

4 

lnP(xj0) = lnP(Y‘^’‘‘’‘|Y(0))-H^lnP(D)'’‘‘’‘|Dj(6>)). (7) 

4=0 
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The dataset is assigned statistical uncertainties 

that are assumed to be normally distributed and uncorre¬ 
lated. Thus, the likelihood lnP(Y'^’^*®'|Y(0)) is modeled 
as a multivariate Gaussian with diagonal covariance. The 
systematic uncertainties on the dataset and their 

effect on the unfolding result are discussed in Sec. mi 

The DC At data sets, in contrast, each comprise a his- 
togrammed distribution of integer-valued entries, and the 
likelihood lnP(D^®'‘^|Dj(0)) is thus more appropriately 
described by a multivariate Poisson distribution. How¬ 
ever, the likelihood calculation for the DCAt data sets 
requires three additional considerations. First, there are 
significant background contributions from a variety of 
sources, as discussed in Section |III E| Secondly, detec¬ 
tor acceptance and efficiency effects are not explicitly ac¬ 
counted for in the DCAt distributions. This implies that 
the total measured yield of signal electrons in each DCAt 
histogram is below what was actually produced, and con¬ 
sequently the measured distributions do not match 

the predictions in normalization. Lastly, because of the 
high number of counts in the region near DCAt = 0, this 
region will dominate the likelihood and be very sensitive 
to systematic uncertainties in the DCAy shape there, 
even though the main source of discrimination between 
charm and bottom electrons is at larger DCAy. 

To deal with the first issue, the relatively normalized 
background described in Sec. |IIIE] is added to each pre¬ 
diction of the DCAt distribution for summed electrons 
from charm and bottom hadrons so that the shape and 
relative normalization of the background component of 
the measurement is accounted for. 

To handle the second, each prediction plus the back¬ 
ground is scaled to exactly match the normalization of 
Ddata^ In this way, only the shape of the prediction is a 
constraining factor. 

To deal with the third, a 5% uncertainty is added in 
quadrature to the statistical uncertainty when the num¬ 
ber of counts in a given DC At bin is greater than a 
reasonable threshold (which we set at 100 counts). This 
accounts for the systematic uncertainty in the detailed 
DCAt shape by effectively de-weighting the importance 
of the region DCAt ~0 while maintaining the overall 
electron yield normalization (as opposed to removing the 
data entirely). This additional uncertainty also necessi¬ 
tates changing the modeling of lnP(D^'^‘®'|Dj(0)) from a 
Poisson to a Gaussian distribution. We have checked that 
varying both the additional uncertainty and the thresh¬ 
old at which it is added has little effect on the results. 


are used to create a decay matrix between hadron pT (ptj 
representing charm hadron px, Pt, or bottom hadron pT, 
Pt) and electron px (Px) and DCAt- Here we treat the 
feed down decay P —>■ P —?> e as a bottom hadron decay 
and exclude it from charm hadron decays. 



[GeV/c] 


4- Decay model and matrix normalization 

The pythia-6 m generator with heavy flavor pro¬ 
duction process included, via the parameter MSEL=4(5), 
is used to generate parent charm (bottom) hadrons and 
their decays to electrons. Electrons within \ri\ < 0.35 
decayed from the ground state charm hadrons (P^, P*^, 
Ds, and Ac) or bottom hadrons P°, Pg, and A;,) 


FIG. 9. (Color Online) (a) The decay matrix, encod¬ 

ing the probability for charmed hadrons decaying to electrons 
within IpI < 0.35 as a function of both electron px (Pt) and 
charm hadron px (px)- (b) An example decay matrix, 
encoding the probability for charmed hadrons decaying to 
electrons within It;! < 0.35 and 1.5 < Px [ GeV/c] < 2.0 as a 
function of both electron DCAt and charm hadron px (Px)- 
In both cases the color intensity represents the probability of 
decay in the given bin. 
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The probability for a charm or bottom hadron at a 
given py to decay to an electron at a given and DCAt 
is encoded in the multidimensional matrices and 

An example decay matrix for charmed hadrons 


is shown in Fig. [9 
spond to the same 


Note that the 17 bins in corre- 
Dins shown along the x-axis in Fig .|15[ 
and that the binning in p^ and DCAt seen in Fig. P is 
the same as that shown in Fig. and Fig. respec¬ 
tively. Furthermore, note that the marginal probabili¬ 
ties do not integrate to unity in these matrices. This 
is because the decay probabilities are normalized to the 
number of hadrons that are generated at all momenta, in 
all directions, and over all decay channels. The probabil¬ 
ity distribution for a hadron integrated over all rapidities 
and decay channels within a given pf^ range to decay to 
an electron at \y\ < 0.35 with a given p^ (integrated over 


DCA-r) is shown in Fig. 10 for an example set of bins. 



FIG. 10. (Color Online) The probability for (a) charm and 
(b) bottom hadrons in a given range of hadron px {jPt and 
Py for charm and bottom hadrons respectively) to decay to 
electrons at midrapidity as a function of electron px (px)- 


hadrons are unmodihed. BRAHMS found that the pion 
and proton Raa did not depend strongly on rapidity up 
to y ~ 3 [SS], justifying the assumption. This assumption 
will further lead us to quote charm and bottom hadron 
yields as a function of pr integrated over all rapidity. The 
second assumption is that all ground state charm hadrons 
experience the same modification as a function of p^- 
While different than the charm suppression, all bottom 
hadrons are assumed to experience the same modifica¬ 
tion. 


An enhancement in the baryon to meson production 
ratios in both nonstrange and strange hadrons has been 
measured at RHIC |59j . which may carry over into the 
heavy quark sector, invalidating the second assumption. 
While there are some models [60j that attempt to in¬ 
corporate this anomalous enhancement into the charm 
hadrons to help explain the measured heavy flavor elec¬ 
tron Raa, there are few measurements to help con¬ 
strain this proposed enhancement. Following Ref. |61) . 
we have tested the effect of this assumption by apply¬ 
ing the observed baryon/meson enhancement to both the 
Kc/D and Ai,/B ratios. As in Ref. [ST], we assume that 
the modification asymptotically approaches 1 for hadron 
Pt > 8 GeV/c. We find that including the enhance¬ 
ment gives a lower charm hadron yield at high-pr and a 
larger bottom hadron yield at high-p^, but the modifi¬ 
cations are within the systematic uncertainties discussed 
in Sec. ITTTh] and shown in Fig. We also find a larger 
bottom electron fraction, which is again within the sys¬ 
tematic uncertainties shown in Fig. 17 While we have 
not used other particle generators to create alternate de¬ 
cay matrices, we find that the and meson px and 
rapidity distributions from pythia are similar to those 
given by Fixed Order -|- Next-to-Leading Log (fonll) 
calculations [33]. We have not included any systematic 
uncertainty due to this model dependence in the final 
result. 


5. Regularization/prior 


To penalize discontinuities in the unfolded distribu¬ 
tions of charm and bottom hadrons, we include a regu¬ 
larization term to the right hand side of equation In 
this analysis we included a squared-exponential function 


In principle, this decay matrix introduces a model de¬ 
pendence to the result. In the creation of the decay ma¬ 
trix we are integrating over all hadron rapidities as well 
as combining a number of hadron species and their de¬ 
cay kinematics to electrons. This involves two assump¬ 
tions. The first is that the rapidity distributions of the 


ln7r(0) = -a^ (ILR^p -f \LR,\^) (8) 


where Rc and Rf, are ratios of the charm and bottom 
components of the parent hadron px vector to the corre¬ 
sponding 17 components of the prior, ©prior, and L is a 
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17-by-17 second-order finite-difference matrix of the form 


/-I 1 \ 

1 -2 1 
1 -2 1 



1 -2 1 
1 -2 1 

\ 1 -1 / 


Thus the addition of this term encodes the assumption 
that departures from ©prior should be smooth by penaliz¬ 
ing total curvature as measured by the second derivative. 

Here, a is a regularization parameter set to a = 1.0 in 
this analysis. We determine a by repeating the unfold¬ 
ing procedure, scanning over a and choosing the value 
of a which maximizes the resulting sum of Eq. and 
— (|LRcp -I- |LRf,p) (Eq. 1^ dropping a^). In this way 
we can directly compare log likelihood values for unfold¬ 
ing results with different a values. We include variations 
on a in the systematic uncertainty as described in Sec¬ 
tion Emi 

We set ©prior to PYTHIA charm and bottom hadron 
Pt distributions scaled by a modified blast wave calcu¬ 
lation [2S] which asymptotically approaches Raa values 
of 0.2(0.3) for D{B) mesons at high-pT- We have tested 
the sensitivity of the result to ©prior by alternatively using 
unmodified pythia charm and bottom hadron px distri¬ 
butions. We find that the result is sensitive to the choice 
of ©prior dominantly in the lowest charm hadron px bins, 
where there is minimal constraint from the data. We have 
included this sensitivity in the systematic uncertainty as 
discussed in Section UlI HI 


6. Parent charm and bottom hadron yield and their 
statistical uncertainty 

The outcome of the sampling process is a distribution 
of © vectors, which is 34-dimensional in this case. In prin¬ 
ciple, the distribution of © vectors contains the full proba¬ 
bility, including correlations between the different param¬ 
eters. The 2-D correlations are shown in Eig. El While 
it is difficult to distinguish fine details in the 34 x 34- 
dimensional grid of correlation plots, we can see a few 
gross features. A circular contour in the 2-D panels rep¬ 
resents no correlation between the corresponding hadron 
Px bins. An oval shape with a positive slope indicates 
a positive correlation between corresponding bins, and 
an oval shape with a negative slope represents an anti¬ 
correlation between corresponding bins. A large positive 
correlation is seen for adjacent bins for high-p^ charm 
hadrons and low-p^ bottom hadrons. This is a conse¬ 
quence of the regularization, which requires a smooth px 
distribution, and is stronger at the higher and lower px 


regions where there is less constraint from the data. We 
also see that, while there is little correlation between the 
majority of nonadjacent px bins, there does seem to be a 
region of negative correlation between the mid to high px 
charm hadrons and the low to mid px bottom hadrons. 
Charm and bottom hadrons in these regions contribute 
decay electrons in the same px region, and appear to 
compensate for each other to some extent. An exam¬ 
ple of this is shown between 3.5 < Px GeV/c < 4.0 and 
2.5 < p^ GeV/c < 3.0 in Fig. [^b)-(d). 

To summarize p(©|x), we take the mean of the 
marginalized posterior distributions (the diagonal plots 
in Fig. 111 for each hadron px bin as the most likely 
values, and the 16*^ and 84*^ quantiles to represent the 
±1 ct uncertainty in those values due to the statistical 
uncertainty in the data modified by the regularization 
constraint. 


7. Re-folded comparisons to data 


The vector of most likely hadron yields, with uncer¬ 
tainties, can be multiplied by the decay matrix to check 
the consistency of the result with the measured data 
(here referred to as re-folding). Figure]^ shows the mea¬ 
sured heavy flavor electron invariant yield in Au-l-Au col¬ 
lisions [12] compared with the re-folded electron spectra 
from charm and bottom hadrons. We find good agree¬ 
ment between the measured data and the electron spec¬ 
trum from the re-folded charm and bottom hadron yields. 
Figure [l^ shows the comparison in electron DGA^ space 
for each bin in electron px- Shown in each panel is the 
measured DGA^ distribution for electrons, the sum of 


the background contributions discussed in Section HIE 


the DCAr distribution of electrons from charm hadron 
decays, and the DGAy distribution of electrons from bot¬ 
tom hadron decays. Note that the sum of the background 
contributions is fixed in the unfolding procedure, and 
only the relative contribution of charm and bottom elec¬ 
trons within IDGAtI < 0.1 cm, as well as their DCAr 
shape, vary. For convenience, the region of the DCAt 
distribution considered in the unfolding procedure is also 
shown, as discussed in Section |HID 6| The sum of the 
background contributions, charm, and bottom electrons 
is shown for a direct comparison with the data. 

The summed log likelihood values for each of the 
DCAt distributions and the electron invariant yield are 
given in Table |lT| To aid in the interpretation of the like¬ 
lihood values, we use a Monte-Carlo method to calcu¬ 
late the expected likelihood from statistical fluctuations 
around the re-folded result. We draw samples from the 
re-folded result based on the data statistics and calcu¬ 
late the distribution of resulting likelihood values. The 
number of standard deviations from the expected value 
is also shown in Table [Hj We find that the log likelihood 
values are large compared to expectations in the heavy 
flavor electron invariant yield as well as the lowest two 
DCAt Px bins. We note that the likelihood values do 
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FIG. 11. (Color Online) The joint probability distributions for the vector of hadron yields, G, showing the 2-D correlations 
between parameters. The diagonal plots show the marginalized probability distributions for each hadron pt bin (i.e. the 
1-dimensional projection over all other parameters). Along the Y-axis the plots are organized from top to bottom as the 17 
charm hadron pr (Pt) bins from low to high p^ followed by the 17 bottom hadron pr (pr) bi ns f rom low to high py- The 


X-axis is organized similarly from left to right. The p^ and Pt binning follows that shown in Fig. 15 The region of green plots 


(top left quadrant) shows the charm hadron yields and the correlations between charm hadron yields. The region of blue plots 
(bottom right quadrant) shows the bottom hadron yields and correlations between bottom hadron yields. The region of orange 
plots (bottom left quadrant) shows the correlations between charm and bottom hadron yields. Sub-panels (b)-(d) show a set 
of example distributions, (b) The 1-D probability distribution of charm hadron yield in 3.5 < py GeV/c < 4.0. (d) The 1-D 
probability distribution of bottom hadron yield in 2.5 < py GeV/c < 3.0. (c) The correlation between (b) and (d). 


not incorporate the systematic uncertainties on the data, In particular the statistical uncertainties on the heavy 
which are handled separately as described in Sec. mE flavor electron invariant yield are much smaller than the 
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systematics at low-pT, making the likelihood value not 
surprising. We find reasonable agreement within uncer¬ 
tainties between the remaining DCAt^ pt bins. 



[GeV/c] 

FIG. 12. (Color Online) The heavy flavor electron invariant 
yield as a function of pr from measured data m compared to 
electrons from the re-folded charm and bottom hadron yields. 
The boxes represent the point-to-point correlated uncertain¬ 
ties on the measured heavy flavor electron invariant yield, 
while the error bars on the points represent the point-to- 
point uncorrelated uncertainties. The label “PHENIX Run 
4 -I- Run 11” on this and all subsequent plots indicates that 
the unfolding result uses the heavy flavor electron invariant 
yield as a function of pr from data taken in 2004 (Run 4) 
combined with DCAt measurements from data taken in 2011 
(Run 11). 


H. Systematic uncertainties 

When performing the unfolding procedure, only the 
statistical uncertainties on the electron DCAj’ and pp 
spectra are included. In this section we describe how we 
consider the systematic uncertainties on both the mea¬ 
sured data and the unfolding procedure. We take the 
following uncertainties into account as uncorrelated un¬ 
certainties: 


TABLE II. The log likelihood values (LL) summed over each 
DCAt distribution and for the comparison to the heavy flavor 
electron invariant yield. Also quoted is the number of data 
points (Np) and the deviation from the log likelihood value 
expected from statistical fluctuations (ALL), as discussed in 
the text, for each comparison. 


Data set 

Np 

LL 

ALL [cr] 

e DCAt 1.5 < Pr < 2.0 

50 

-195.5 

-3.8 

e DCAt 2.0 < Pp < 2.5 

50 

-156.5 

-2.9 

e DCAt 2.5 < pp < 3.0 

50 

-115.8 

-0.6 

e DCAt 3.0 < pp < 4.0 

50 

-104.1 

-1.8 

e DCAt 4.0 < pp < 5.0 

50 

-53.2 

0.0 

e Inv. Yield. 1.0 < Pt < 9-0 

21 

-45.9 

-3.5 

Total Sum 

271 

-673.8 



1. Systematic uncertainty in the heavy flavor electron 
Pp invariant yield 

2. Uncertainty in the high-multiplicity background 

3. Uncertainty in the fraction of nonphotonic elec¬ 
trons (Lnp) 

4. Uncertainty in normalization 

5. Regularization hyperparameter a 

6. Uncertainty in the form of 0prior 


The uncertainty in Fnp (See Sec. VIA), and ATeS are 
propagated to the unfolded hadron yields by varying each 
independently by Tier, and performing the unfolding pro¬ 
cedure with the modified background template. The dif¬ 
ference between the resulting hadron yields and the cen¬ 
tral values is taken as the systematic uncertainty. The 
same procedure is used to determine the uncertainty in 
the result due to the regularization parameter, which is 
varied by (’(q ™ based on where the summed likelihood 
from both the data and regularization drops by 1 from 
the maximum value. 

The uncertainty in the high-multiplicity background 
includes two components. The first is the uncertainty 
on the normalization of the high-multiplicity background 
DCAj- distribution, as determined in Sec. |III E 2| and 
shown in Fig. This is propagated to the unfolded 
hadron yields by varying the normalization by ±lcr and 
performing the unfolding procedure with the modified 
background template, as with the Fnp and ATeS uncer¬ 
tainties. The second component addresses the small ex¬ 
cess in the embedded primary electron distribution ob¬ 
served in Fig.j^and not accounted for by using the DCAt 
distribution for large DCAp,. We parametrize the excess, 
which is more than two orders of magnitude below the 
peak, and apply it to the background components, re¬ 
performing the unfolding procedure to find its effect on 
the hadron yield. Both effects combined are small rela¬ 
tive to the dominant uncertainties. 
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Unfolded bottom 
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FIG. 13. (Color Online) The DCAj distribution for measured electrons compared to the decomposed DCAj distributions 
for background components, electrons from charm decays, and electrons from bottom decays. The sum of the background 
components, electrons from charm and bottom decays is shown as the red (upper) curve for direct comparison to the data. The 
gray band indicates the region in DCAj considered in the unfolding procedure. Also quoted in the figure is the bottom electron 
fraction for |DCAt| <0.1 cm integrated over the given pr range. The legend follows the same order from top to bottom as 
panel (b) at DCAj = —0.1 cm. 


Incorporating the pt correlated systematic uncertainty ficult. Ideally one would include a full covariance ma- 
on the heavy flavor electron invariant yield is more dif- 
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trix encoding the pr correlations into the unfolding pro¬ 
cedure. In practice, the methodology employed in [12] 
does not provide a convenient description of the pt cor¬ 
relations needed to shape the covariance matrix. Instead 
we take a conservative approach by considering the cases 
which we believe represent the maximum px correlations. 
We modify the heavy flavor electron invariant yield by ei¬ 
ther tilting or kinking the spectrum about a given point. 
Tilting simply pivots the spectra about the given point 
so that, for instance, the first point goes up by a fraction 
of the systematic uncertainty while the last point goes 
down by the same fraction of its systematic uncertainty, 
with a linear interpolation in between. Kinking simply 
folds the spectra about the given point so that that the 
spectrum is deformed in the form of a V. We implement 
the following modifications and re-perform the unfolding 
procedure: 

1. Tilt the spectra about px = 1.8 GeV/c by Tier of 
the systematic uncertainty. 

2. Tilt the spectra about px = 5 GeV/c by ±la of 
the systematic uncertainty. 

3. Kink the spectra about px = 1.8 GeV/c by Tier of 
the systematic uncertainty. 

4. Kink the spectra about px = 5 GeV/c by iltr of 
the systematic uncertainty. 
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The Px points about which the spectra were modihed 
were motivated by the points in px at which analysis 
methods and details changed, as discussed in [12]. We 
then take the RMS of the resulting deviations on the 
hadron yield from the central value as the propagated 
systematic uncertainty due to the systematic uncertainty 
on the heavy flavor electron invariant yield. 

The effect of our choice of 0prior on the charm and 
bottom hadron yields is taken into account by varying 
^priori as discusscd in Section |III G 5] The differences 


between each case and the central value are added in 
quadrature to account for the bias introduced by ©prior- 

The uncertainties on the unfolded hadron yields due to 
the six components described above and the uncertainty 
determined from the posterior probability distributions 
are added in quadrature to give the uncertainty shown in 

Fig.jTsj 

Due to the correlations between charm and bottom 
yields, the relative contributions from the different un¬ 
certainties depend on the variable being plotted. To give 
some intuition for this, we have plotted the relative con¬ 
tributions from the different uncertainties to the fraction 
of electrons from bottom hadron decays as a function of 
Px (discussed in Sec. IV A) in Fig. 14 One can see that 
the dominant uncertainties come from the statistical un¬ 
certainty on the DCA-r and heavy flavor electron invari¬ 
ant yield, the systematic uncertainty on the heavy flavor 
electron invariant yield, and Fnp- We remind the reader 


FIG. 14. (Color Online) The relative contributions from the 
different components to the uncertainty on the fraction of 
electrons from bottom hadron decays as a function of px- 
The shaded red band in each panel is the total uncertainty. 


that for Px > 5 GeV/c we no longer have DCAt- infor¬ 
mation to directly constrain the unfolding, and all infor¬ 
mation comes dominantly from the heavy flavor electron 
invariant yield, leading to the growth in the uncertainty 
band in this region. 


IV. RESULTS 


The final result of the unfolding procedure applied si¬ 
multaneously to the heavy flavor electron invariant yield 
vs Px (shown in Fig. 12) and the hve electron DGA-p 
distributions (shown in Fig. 13) is the invariant yield of 


charm and bottom hadrons, integrated over all rapidity, 
as a function of px- As a reminder, the hadron yields 
are integrated over all rapidity by assuming the rapid¬ 
ity distribution within pythia is accurate and that it is 
unmodified in Au-|-Au, as detailed in Sec. jlll G 4| The 
unfolded results for MB (0%-96%) Au-|-Au collisions at 
=200 GeV are shown in Fig. 
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The central point 


represents the most likely value aiid the shaded band 
represents the Icr limits on the combination of the un- 
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certainty in the unfolding procedure and the systematic 
uncertainties on the data, as described in Sec. Ill H The 
uncertainty band represents point-to-point correlated un¬ 
certainties, typically termed Type B in PHENIX publi¬ 
cations. There are no point-to-point uncorrelated (Type 
A), or global scale uncertainties (Type C), from this pro¬ 
cedure. 


The uncertainties on the hadron invariant yields shown 
in Fig. grow rapidly for charm and bottom hadrons 
with pt > 6 GeV/c. This is due to the lack of DCAy in¬ 
formation for pf, > 5 GeV/c. Above pf. > 5 GeV/c, the 
unfolding is constrained by the heavy flavor electron in¬ 
variant yield only. This provides an important constraint 
on the shape of the hadron pr distributions, but the 
DCAt distributions provide the dominant source of dis¬ 
criminating power between the charm and bottom. How¬ 
ever, due to the decay kinematics, even high pT hadrons 
contribute electrons in the range 1.5 < pf. [ GeV/c] < 
5.0. We find that charm(bottom) hadrons in the range 
7 < p^ [ GeV/c] < 20 contribute 18.2%(0.3%) of the to¬ 
tal electron yield in the region 1.5 < pf. [ GeV/c] < 5.0. 
This explains the larger uncertainties in the bottom 
hadron yield compared to the charm hadron yield at high 
Pt- 

The yield of mesons over |p| < 1 as a function of 
Pt has been previously published in Au-|-Au collisions 
=200 GeV by STAR [T3]. In order to com¬ 
pare our unfolded charm hadron results over all rapidity 
to the STAR measurement, we use pythia to calculate 
the fraction of mesons within \y\ < 1 compared to 
charm hadrons over all rapidity. Since the measurement 
by STAR is over a narrower centrality region (0%-80% 
vs 0%-96%), we scale the STAR result by the ratio of the 
A^coii values. This comparison is shown in Fig. |16| For 
added clarity, we have fit the STAR measurement with a 
Levy function modified by a blast wave calculation given 
by 


^ _ (1 -pi)pt \ 
P2 ) 




Upt) = Po 

X f 1.3^27rp^G(pr,P3,P4) + 


( 10 ) 


P5 


1 -I- 


where G(jpt,P^,Pa) is a standard Gaussian function, and 
Pi are the parameters of the fit. The ratio of the data to 
the fit is shown in the bottom panel of Fig. We find 
that, within uncertainties, the unfolded 11° yield agrees 
with that measured by STAR over the complementary 
Pt range. The unfolded yield hints at a different trend 
than the STAR data for p^ > 5 GeV/c. However, we 
note that the {pt) of charm(bottom) hadrons which con¬ 
tribute electrons in the range 4.0 < pt [ GeV/c] < 5.0 
is 7.2(6.4) GeV/c. This means that the yields of charm 
and bottom hadrons have minimal constraint from the 
DCAt’ measurements in the high-pT’ regions, which is 
represented by an increase in the uncertainties. 


A. The bottom electron fraction 


The fraction of heavy flavor electrons from bottom 
hadrons () is computed by re-folding the charm 
and bottom hadron yields shown in Fig. to get the in¬ 
variant yield of electrons from charm and bottom decays 
at midrapidity (|p| < 0.35). Here the electrons from bot¬ 
tom hadron decays include the cascade decay 6 —?► c —>■ e. 
The resulting bottom electron fraction is shown as a func¬ 


tion of Pt in Fig. 17 The central values integrated over 
the Pt range of each DCA-r distribution are also quoted 
in Fig. As in the hadron yields, the band represents 
the Icr limits of the point-to-point correlated (Type B) 
uncertainties. 

Also shown in Fig. [l^is the bottom electron fraction 
predictions from fonll [33] for p-l-p collisions at 
=200 GeV. We find a bottom electron fraction which 
is encompassed by the fonll calculation uncertainties. 
The shape of the resulting bottom electron fraction shows 
a steeper rise in the region 2.0 < pt [GeV/c] < 4.0 with a 
possible peak in the distribution compared to the central 
FONLL calculation. 

The fraction of electrons from bottom decays has 
been previously measured in p-l-p collisions at y/s~^ 
=200 GeV by both PHENIX [31] and STAR [^. 
These measurements are made through electron-hadron 
or electron-D meson correlations. These are very differ¬ 
ent analyses than the one presented here, and have their 


own model dependencies. In Fig. 18 we compare the 
bottom electron fraction between our unfolded Au-|-Au 
result and the electron-hadron correlation measurements 
in p+p. For p-r > 4 GeV/c we find agreement between 
Au-|-Au and p-l-p within the large uncertainties on both 
measurements. This implies that electrons from bottom 
hadron decays are similarly suppressed to those from 
charm. For reference, included in Fig. [^ is the central 
FONLL calculation which, within the large uncertainties, 
is consistent with the p+p measurements. 

With the additional constraints on the bottom elec¬ 
tron fraction in p+p from the correlation measurements 
and the measured nuclear modification of heavy flavor 
electrons, we can calculate the nuclear modification of 
electrons from charm and bottom hadron decays sepa¬ 
rately. The nuclear modifications, R'Xa^ 
charm and bottom hadron decays respectively are calcu¬ 
lated using 


R 

R 


c—¥e _ 

AA — 

b—>-e 

AA 


(1-FauAu) dHF 


FauAi 

p 

JTpp 


R 


HF 

AA) 


( 11 ) 

( 12 ) 


where FauAu and Fpp are the fractions of heavy flavor 
electrons from bottom hadron decays in Au-|-Au and p-l-p 
respectively and R^^ 1® nuclear modification of heavy 
flavor electrons (combined charm and bottom). Rather 
than combining all measurements for the bottom electron 
fraction in p+p, which introduces a further extraction un¬ 
certainty, we have chosen to calculate R'^^ and R!’^ us- 
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FIG. 15. (Color Online) Unfolded (a) charm and (b) bottom hadron invariant yield as a function of pr, integrated over all 
rapidities, as constrained by electron yield vs DCAt in 5 Px bins and previously published heavy flavor electron invariant yield 
vs px |12| . 


ing only the six STAR electron-hadron Fpp values. When 
performing the calculation we determine the full proba¬ 
bility distributions assuming Gaussian uncertainties on 
.PauAu, Fpp and when determining the charm 

and bottom hadron yields, we take the median of the dis¬ 
tribution as the central value, and the 16% and 84% of 
the distribution as the lower and upper la uncertainties. 
The resulting values are shown in Fig. l^a). We find 
that the electrons from bottom hadron decays are less 
suppressed than electrons from charm hadron decays for 
3 < pt GeV/c < 4. To further clarify this statement we 
calculate the ratio of Raa^/F-aT^ shown in Fig. [T^b). 
In this ratio, the uncertainty on cancels. Here again 
we calculate the full probability distributions and use the 
same procedure as above to determine the central values 
and uncertainties. We find that the probability distri¬ 
butions for R^aa^/R'^^ are highly nonGaussian, which 
leads to the large asymmetric uncertainty band shown 
in Fig. &). It is clear from the ratio that 6 —>■ e is 
less suppressed than c —?> e at the la level up to pt ~ 4 
GeV/c. 


V. DISCUSSION 

There are a number of theoretical calculations in the 
literature for the interaction of charm and bottom quarks 
with the QGP. Many of these models have predictions for 


the nuclear modification factor Raa for electrons from 
charm decays and, separately, Raa for electrons from 
bottom decays. For consistency, we have assumed the 
FONLL [33] yields for electrons from charm (bottom) de¬ 
cays calculated for p+p at ^/s^ =200 GeV and then 
scaled them by the heavy-ion model results for the Raa 
of electrons from charm (bottom). 

Figure 20 a) compares the bottom electron fraction 
from one class of calculations modeling only energy loss 
of these heavy quarks in medium. In an early pQGD cal¬ 
culation by Djordjevic, Gyulassy, Vogt, and Wicks [33], 
the authors apply the DGLV theory of radiative energy 
loss. They find that even for extreme opacities with gluon 
rapidity densities up to 3500, the bottom quark decay 
electrons dominate at high-pj’ and that limits the single 
electron Raa to the range 0.5-0.6 for px > 5 GeV/c. 
Although this result is known to be higher than the 
PHENIX measured heavy flavor electron Raa [E], we 
show the b -A e/{b -A e + c ^ e) predictions for gluon 
rapidity densities of 1000 and 3500 in Fig. 20'a). How¬ 
ever, we do note that the calculations are for 0%-10% 
central collisions compared to the MB data, although 
the calculations span a factor of 3.5 range in the gluon 
density. We find that the calculations for both gluon ra¬ 
pidity densities are in good agreement with our results 
for Pt < 4 GeV/c, but are slightly above and outside 
the uncertainty band on the unfolded result at higher 
Pt- More recent calculations in the same framework, but 
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FIG. 16. (Color Online) The invariant yield of mesons 
as a function of pr for \y\ < 1 inferred from the nnfolded 
yield of charm hadrons integrated over all rapidity compared 
to measurements from STAR m- See the text for details 
on the calculation of the yield inferred from the unfolded 
result. To match the centrality intervals, the STAR result has 
been scaled by the ratio of Acoii values. The bottom panel 
shows the ratio of the data to a fit of the STAR yield. 



FIG. 17. (Color Online) The fraction of heavy flavor electrons 
from bottom hadron decays as a function of pt from this work 
and from fonll p+p calculations m - 



with the inclusion of collisional energy loss [31] , result in 
a heavy flavor electron high-p^ Raa closer to 0.3 and in 
reasonable agreement with previous PHENIX published 
results [12] . This updated prediction for the bottom elec¬ 
tron fraction, also shown in Fig. gives a similar value 
to their previous result, but is only published for > 5 
GeV/c. 

Figure [^b) compares the bottom electron fraction 
from a calculation using a T-matrix approach by van 
Hees, Mannarelli, Greco, and Rapp [53] . The authors 
provided us with different results for 0%-10% central 
Au-|-Au collisions depending on the coupling of the 
heavy-quark to the medium. The coupling is encapsu¬ 
lated in the diffusion parameter D, where smaller values 
yield a stronger coupling. Shown in Fig. [^b) are three 
results corresponding to three values of the parameter 
D{2ttT) = 4,6,30. The largest D value, correspond¬ 
ing to the weakest coupling, yields almost no deviation 


FIG. 18. (Color Online) bottom electron fraction as a func¬ 
tion of Pt compared to measurements in p+p collisions at 
=200 GeV from PHENIX [34] and STAR [35]. Also shown are 
the central values for fonll [33] for p+p collisions at 
=200 GeV. 


from the p+p reference fonll result, and the successively 
stronger coupling pushes the bottom fraction contribu¬ 
tion higher and higher. We find that the calculations 
with D{2nT) =4,6 are in good agreement with our re¬ 
sult for Pt < 4 GeV/c, but begin to diverge where the 
calculation stops at 5 GeV/c. 

Figure 20 c) compares the bottom electron fraction 
from another class of calculations which employ a combi¬ 
nation of Langevin, or transport type modeling of heavy- 
quarks, in the bulk QGP with energy loss mechanisms 
that dominate at higher pp- In Ref. [64], Alberico et al. 
employ a Langevin calculation where a good match to 
the PHENIX heavy flavor electrons is found. It is no- 
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[GeV/c] 


FIG. 19. (Color Online) (a) The Raa for c ^ e, b ^ e and 
combined heavy flavor m as a function of ■ The c —>■ e and 
6 —>■ e Raa is calculated using Eq.jlTJ^ where FauAu uses the 
unfolded result determined in this work and Fpp determined 
from STAR e — h correlations [35]. (b) The ratio R!aa^/R’aa^ 
as a function of p^. 


table that this calculation has a very strong suppression 
of charm decay electrons such that bottom contributions 
dominate even at modest pt > 2 GeV/c. The calcu¬ 
lations are consistent with the data for p-r < 4 GeV/c 
and over-predict the bottom contribution for higher px 
values. 

Figure 20 c) also compares the bottom electron frac¬ 
tion from another variant of the Langevin calculation by 
Gao et al, as detailed in Ref. [HS]. For this calculation, 
we show two results corresponding to two different input 
values D{2'kT) = 1.5 and 6. For the lower parameter, 
again stronger heavy-quark to medium coupling, there 
is a sharp rise in the bottom contribution which then 
flattens out. This feature is due to the increased colli- 
sional energy loss, which has a larger effect on the charm 
quarks, coupled with the strong radial flow effects en¬ 
abling the heavier bottom quarks to dominate even at 
Pt ^ 2 GeV/c. These calculations use an impact pa¬ 
rameter of 6 = 6.5 fm, which should roughly correspond 
to MB collisions. We find that the calculation using the 


larger value of D{2'kT) = 6.0 is in reasonable agreement 
with the data across the calculated px range. 

Lastly, Fig. I^d) shows a more recent calculation by 
He et al. employing a T-matrix approach similar to that 
shown in Fig. [20](b), but with a number of updates as 
described in Ref. |66|. In this case the authors provided 
a calculation of the bottom electron fraction in both p-l-p 
and Au-I-Au at =200 GeV, and we therefore do not 

calculate the bottom fraction using fonll as a baseline. 
The calculation is performed for the 20%-40% centrality 
bin, which the authors find well represents MB. We find 
that the calculation under-predicts the bottom fraction 
for pt < 3 GeV/c, although it is worth noting that the 
calculation in p-|-p is also below the fonll curve across 
the full px range. Above p^ ~ 3 GeV/c the calculation 
is in agreement with the measurement. It is also worth 
noting that, of the models presented here, this is the only 
one that shows in Au-I-Au a slight decrease in the bottom 
fraction at high px ■ 

There are numerous other calculations in the litera¬ 
ture [67H6^ that require mapping charm and bottom 
hadrons to electrons at midrapidity to make direct data 
comparisons. We look forward to soon being able to test 
these calculations with analysis of new PHENIX data 
sets. 


VI. SUMMARY AND CONCLUSIONS 

This article has detailed the measurements of electrons 
as a function of DGAt and px from Au-I-Au data taken 
=200 GeV in 2011 with the enhanced vertexing 
capabilities provided by the VTX detector. In conjunc¬ 
tion with previous PHENIX results for the heavy flavor 
electron invariant yield as a function of px m , we per¬ 
form an unfolding procedure to infer the parent charm 
and bottom hadron yields as a function of px- We find 
that this procedure yields consistent agreement between 
the heavy flavor electron invariant yield and the newly 
measured electron DCAt distributions. 

We find that the extracted 19° yield vs px is in good 
agreement with that measured by STAR m over the 
complimentary px region. Without a proper p-l-p base¬ 
line extracted from a similar analysis it is difficult to 
make any quantitative statements about the charm or 
bottom hadron modification. 

We compare the extracted bottom electron fraction to 
measurements in p-l-p collisions and find agreement be¬ 
tween Au-I-Au and p+p for px > 4 GeV/c within the 
large uncertainties on both measurements. The agree¬ 
ment between Au-I-Au and p-fp coupled with the mea¬ 
sured heavy flavor electron Raa strongly implies that 
electrons from charm and bottom hadron decays are 
suppressed. Using these components we calculate the 
nuclear modification for electrons from charm and bot¬ 
tom hadron decays and find that electrons from bottom 
hadron decays are less suppressed than those from charm 
hadron decays in the range 3 < px GeV/c < 4. We fur- 
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FIG. 20. (Color Online) Bottom electron fraction as a function of pr compared to a series of model predictions detailed in the 
text. 


ther compare the bottom electron fraction to a variety 
of model calculations employing variously energy loss, 
Langevin transport, and T-matrix approaches. We find 
that there are a number of models which are in reasonable 
agreement with the extracted bottom electron fraction 
within the relatively large uncertainties. 

We note that a significantly larger data set of Au+Au 
collisions at =200 GeV was collected in 2014 with 

an improved performance of the VTX detector. The 2014 
Au+Au data coupled with the p+p data taken in 2015 
should yield both an important baseline measurement of 
the bottom electron fraction and a more precise measure¬ 
ment in Au+Au. 
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APPENDIX: DETAILED NORMALIZATION OF 
ELECTRON BACKGROUND COMPONENTS 

This appendix details the calculation of the normaliza¬ 
tions for the background components: 

• Photonic electrons 


Here is the ratio of conversion electrons to elec¬ 

trons from Dalitz decays in the 2004 PHENIX detector 
calculated by a full GEANt3 simulation. The factors 

• rDahtz{Tr°) 

• rDaUtziv) 

• rDahtzil) 

are the fractional contributions of t], and direct pho¬ 
ton contribution to the total Dalitz decays, respectiveljQ 
We only consider the contributions of tt*^, 77 , and jdir (di¬ 
rect photon) since the sum of other contributions is small 
(5% or less). Thus they are normalized such that 


'^rDalitzii) = (15) 


Figure 21 shows rjjautz for 7r°, ij, and direct photon 
as a function of transverse momentum of the electrons 
for MB Au-|-Au collisions at 200 GeV. The ratios are 
calculated from the invariant yield of 7r°[49], 77 ( 50 ], and 
direct photons[5llM]. 


• Kaon decay electrons 

• Heavy quarkonia decay electrons 


using the bootstrap method described in Sec. |HIF| We 
first determine the fraction of nonphotonic electrons, 
Fnp . We then calculate the normalization of Dalitz and 
conversion components followed by the normalization of 
Kes and quarkonia components. 


A. Fraction of nonphotonic electrons Fnp 

We first determine Fnp, the fraction of nonphotonic 
electrons to inclusive electrons after the application of 
all analysis cuts, including the conversion veto cut. Note 
that nonphotonic electrons include contributions from 
heavy flavor semi-leptonic decays, quarkonia decays, and 
kaon decays. Photonic electrons are from 7 r° and 77 Dalitz 
decays and photon conversions. 

Fnp in the 2011 data can be determined using the pub¬ 
lished 2004 result (T^j as follows. Let Ynp be the yield of 
nonphotonic electrons and Yuaiitz the yield of electrons 
from Dalitz decays. Note that both Inp and Yoaiitz are 
independent of the year of data taking. In the PHENIX 
2004 Au-|-Au data run, the ratio of the nonphotonic elec¬ 
tron yield to the photonic electron yield (R^p^) '''^as mea¬ 
sured. The relation of Ynp and Yjjautz is as follows: 

Inp = RnT(1 + Rcd) X YoaUtz. (13) 

where R^^ represents the ratio of conversion electron 
yield to Dalitz electron yield in the 2004 PHENIX detec¬ 
tor. It is calculated as 

RcD= E RcD\^-rDaUtz{^)- (M) 

2=7r®,?7,7 



FIG. 21. (Color Online) The fraction of 7r°, rj, and direct 
photon Dalitz decay electrons in all Dalitz electrons as a func¬ 
tion of electron pr (Pt)- 

In the 2011 data set the observed electron yields from 
conversion and Dalitz decays are modified by the elec¬ 
tron survival probability after the conversion veto cut is 


^ Here we include internal conversion of direct photon in Dalitz 
decays. Note that the Dalitz decay of {rj) is caused by internal 
conversion of one of two decay photons in > 77 . 
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applied. The yield of photonic electrons which pass the 
conversion veto is 

=Rfh^ xYoaUt., (16) 

RfD= E {SD{l) + Sc-R^SDi^)rDamS){^^) 

,77,7 

where Sc is the survival probability of conversion elec¬ 
trons, Soiri), Soil) are survival probabilities of 

Dalitz decay electrons from t], and direct photons, 
respectively, as shown in Fig. (* = 7’'°,??,7)is 

the ratio of conversion electrons to Dalitz electrons for 
particle i in the 2011 PHENIX detector after the addi¬ 
tion of the VTX and the replacement of the beam pipe. 
It is determined to be R^^{i) « 1.10 from full GEANt3 
simulations. 
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FIG. 22. (Color Online) The fraction of nonphotonic elec¬ 
trons to inclusive electrons as a function of electron pr (Pt)- 



p® [GeV/c] 


FIG. 23. (Color Online) The fraction of tt'^, rj, and direct 
photon electrons in all photonic electrons as a function of 
electron Pt (Pt)- 


where is the number of electrons with conversion veto 
after the subtraction of the hadronic contamination and 
random background. 


B. Normalization of Dalitz and conversion 
components 

In the previous section we obtained Np, the number of 
photonic electrons in the data after the conversion veto 
cut. There are two components in the photonic electrons 
(iV|.). 

1. Electrons from Dalitz decays (tt^ - 1 - 77 - 1 - 7 ) 


The fraction of nonphotonic electrons to inclusive elec- 2. Electrons from conversions in the beam pipe and 

trons can then be calculated as ^ 


' Vnp + F|«n- 

_ /iq'i 

i?20p04(l + ^g.04)+^2pn I 

The resulting Fnp as a function of Pp and the calcu¬ 
lated systematic uncertainties due to the uncertainties on 
the input yields is shown in Fig. With Fnp in hand, 
we obtain the number of photonic electrons, Np, and the 
number of nonphotonic electrons, N^p as 


Nf, = Ne{l - Fnp) (20) 

fV^P = NeF^p, (21) 


In the next step, we determine the normalization of 
Dalitz and conversions separately. This is needed since 
the shape of DCA-r distribution of Dalitz and conversions 
are different. 

After application of the conversion veto cut, we have 


iV§(7) = ScRh°hHi)i^ - Srandom)eAYDaUtz, (22) 

~ ^random)^AYDalitz^ (23) 

(f = 7r°,?7,7) (24) 

where N^(i) and Npii) are the number of electrons from 
conversions and Dalitz from particle i after the conver¬ 
sion veto cut, respectively; Srandom is the common reduc¬ 
tion factor of tracks due to random hits in the windows 
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of the conversion veto cut; and tA is the efficiency and 
acceptance without the conversion veto cut. Since the 
number of photonic electron is Np{i) = Np{i) + 
the fraction of conversions and Dalitz decays in the pho- 


tonic electrons are 




ScRcoii) 

(25) 

Nm 

SD{i) + ScR^SD{iy 

Nh{i) 

Snii) 

(26) 

Nm 

SD{i) + ScR^SD{y' 


The fraction of electrons from conversions (N^/Np) 
and Dalitz {Np/Np) is the average of these fractions, 
thus: 


= N’p 


rP^uy _ ScRcDii) 


Nh = N^P E 


Soil) + ScR^^hH^) 
Soji) 

Soit) + ScRl°h\^ 


, (28) 


where rP^{i), {i = 7r°,?7,7) is the relative contribution of 
electrons from (conversion + Dalitz decay) for particle i 
after application of conversion veto cut. Figureshows 
rP^{i) {i = TT'^,ri,'y) as a function of pp. The conversion 
contributions are nearly the same for 7r°,?7 and 7, and 
effectively cancel when calculating the ratio. Therefore, 
pph (Fig. 23) is almost identical with roaiitz (Fig. 21). 


C. Normalization of Ke 3 and quarkonia components 


The ratio of electrons from kaons to all nonphotonic 
electrons before the application of the conversion veto 
cut, Sk, is calculated from the ratio of the nonphotonic 
electron yield to the electron yield from kaons [H] . Com¬ 
pared to Ref. [12], we find that ^ 50% of electrons from 
kaon decays are removed by DC At and DC At cuts as 
well as the method used to subtract random background, 
which contains some real electrons from kaon decays. 


The ratio of electrons from J/'0 decays to all nonpho¬ 
tonic electrons before the application of the conversion 
veto cut, Sj/^, is taken from Ref. |12] . The survival rate 
for electrons from J/ij) decays , Sj/^, is taken to be unity, 
while the survival rate for ATeS decays, Sk, is taken to 
be the same value as that for electrons from charm and 
bottom decays (namely, 5'hf)- See Sec. [Ill E 3|for details. 


After application of conversion veto cut, the normal¬ 
izations of these two nonphotonic electron components 
are described by 


^NP (29) 
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